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ABSTRACT 

We present the results of a variability study of accreting young stellar objects in the Chameleon 
I star-forming region which is based on ~ 300 high resolution optical spectra from the multi- 
object fibre spectrograph FLAMES/GIRAFFE at the ESO/VLT. Twenty five objects with spec- 
tral types from G2 -M5.75 were observed 12 times over the course of 15 months. Using the 
emission lines Ha (6562.81 A) and Call (8662.1 A) as accretion indicators we found 10 ac- 
^ c*j creting and 15 non-accreting objects. 

Q_[ We derived accretion rates for all accretors in the sample using the Ha equivalent width, 

Ha 10% width and the Call (8662.1 A) equivalent width. We found the Ha equivalent width 
of accretors varied by ~ 7 - 100 A over the 15 months. This corresponds to a mean amplitude 
of variations in derived accretion rate of ~ 0.37 dex. The amplitude of variations in derived 
accretion rate from Call equivalent width were ~ 0.83 dex and from Ha 10% width were 
~ 1 . 1 1 dex. Based on the large amplitude of variations in accretion rates derived from the Ha 
10% width with respect to the other diagnostics, we do not consider it to be a reliable accretion 
rate estimator. Taking the variations in Ha EW and Ca II EW accretion rates to be closer to 
the true value, they suggest that the spread which has been found around the accretion rate to 
stellar mass relation is not due to the variability of individual objects on time-scales of weeks 
to ~ 1 year. From these variations we can also infer that the accretion rates are stable within 
< 0.37 dex over time-scales of less than 15 months. 

A major portion of the accretion variability was found to occur on less than the shortest 
time-scales in our observations, 8-25 days, which is comparable with the rotation periods of 
these young stellar objects. This could be an indication that what we are probing is spatial 
structure in the accretion flows, and also suggests that observations on time-scales of ~ a 
couple of weeks are sufficient to limit the total extent of accretion rate variations in typical 
young stars. 

No episodic accretion was observed, all 10 accretors continuously accreted for the entire 
period of observations, and though they may have undetected, low accretion rates, the non- 
accretors never showed any large changes in their emission that would imply a jump in their 
accretion rate. 

Key words: accretion, accretion discs stars: variables: T Tauri - star: formation - stars: pre- 
main-sequence 



1 INTRODUCTION 

Accretion is a vital feature in the process of star formation as it 
controls the flow of matter and angular momentum from the cir- 
cumstellar environment on to the young stellar object (YSO). For 
low-mass YSOs, th e current paradigm for this proc ess is magne- 
tospheric accretion dCamenzin dl ll99d : lK6nigilll99lh : In short, the 
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strong magnetic field of the young star truncates the inner part of 
the disc and channels the acc retion flow from the inner edge of the 
disc to the star (see review bv lBouvier et alj .2007). 

Accretion causes a continuous excess flux from the ultra-violet 
(UV) to the near-infrared as well as a spectrum of emission lines, in- 
cluding the H Balmer series, the Call infrared triplet (Ca IRT), and 
He I lines. Given the relevance of accretion, it is somewhat worri- 
some that our understanding of the observational indicators for ac- 
cretion is poor, in particular their relation to the mass accretion rate 
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M. The analysis of the lines and continuum relies on simplified as- 
sumpti ons about the geometry and structure of the accretion shock 
(e.g. lHerczeg & Hillenbrand 2008). Moreover, many emission lines 
have additional components due to strong winds/outflows, either 
from the star or from the disc dRav et ail2007l) . 

One promising way to improve our understanding is to study 
accretion spectroscopically in the time-domain. The majority of the 
spectroscopic monitoring studies of YSOs have focused on a rela- 
tively small number of targets, specifically selected for their known 
variability. Some of these studies have le d to a detailed close- 
up view on the disc -accretion system (e.g, iGiampapa et alj[l993l : 
lBouvieretaflll999l) and provided unique insig hts into the mag- 
netospheric accretio n/ejection scenario (e.g. Ijohns & Basra Il995l : 
iBouvier et alj 120071) . In other cases the changes in the variabil- 
ity patterns are highly complex and difficult to interpret (e.g. 
lAlencar et aT1l2005h . 

Various models have been put forward to explain accretion 
variability. Misalignment between rotation and stellar magnetic 
field axes will result in a non-uniform accretion flow. Different parts 
of the flow will then come into view as the star rotates, changing the 
emis sion signature we obse rve on the time-scale of the rotation pe- 
riod dKurosawa et al.l l2008). Non-ax isymmetric fields will also in- 
troduce asymmetry into the system dMohantv~& Shu 2008). Addi- 
tionally inhomogeneities in the disc will cause the observed accre- 
tio n rate to vary, as is th e case with an accreting T Tauri star found 
bv lBouvieretalJdl999h with roughly constant brightness, broken 
by a quasi-periodic fading. This was attributed to a wall in the inner 
circumstellar disc which obscures the light from the central source 
periodically as it rotates. In the same way, such a wall could be a 
source of apparent accretion variability on a much longer time-scale 
if plac ed further out in the circumstellar disc. IVorobvov & Basul 
(2010) showed that for young embedded objects, gravitational in- 
stabilities in the disc and disc fragmentation will lead to large scale 
accretion variations, such as FUor and EXor events where the accre- 
tion rates can vary between 10 _4 MQyr _1 and lO _8 M0yr _1 over 
very long periods of thousands of years. Another source of variabil- 
ity is differential rotation between the disc and the star. This can 
lead to expanding field lines, which could radically reduce the mass 
accretion rate as the magnetic field lin es open, and will the n cause 
it to increase again upon reconnection dGoodson et aiTl 998). These 
magnetic reconnection led variations are expected to occur on the 
time -scales on the order of the rotation period of the circumstellar 
disc dHavashietalll996T) . 

Apart from a few benchmark objects, observations have yet to 
distinguish between these scenarios. Each of these cases will show 
different amplitudes and time-scales of variations. In particular vari- 
ability studies can provide infor mation on the spa tial structure and 
dynamics of the accretion flow dBouvier et alj2003h . They can also 
help to understand the origin of the scatter in M at a given stel- 
lar m ass M* , and refine the M ~ hit correlation dMuzerolle et alj 
2003). The scatter within this relation can be up to two orders of 
magnitude for a given mass, and it is not clear if this scatter comes 
from time-varying accretion rates of individual objects, differences 
in the average accretion rates of objects of the same mass or both. 

iNguven et al] (2009) found that the variability in accretion 
rates of a sample of YSOs is dominated by variations on time-scales 
of a few days, and also the amplitude of variations is not enough to 
explain the full spread in the M-M* re lation. This is in agreement 
with previous studies dNatta et alj|2004l) . and suggests that initial 
conditions or environmental differences are more likely to be the 
cause of this spread. Howeve r long term variations could also play 
a role dDullemond et alj2006h . Our observations are designed to in- 



vestigate the variability on longer time-scales than previous studies 
and better define the spread associated with the M-M* relation. 

In this paper, we set out to investigate the extent and char- 
acteristics of accretion-related variability in YSOs in an unbiased 
sample of objects, not chosen for their variability. We observed 
25 low-mass YSOs in the Chamaeleon-I star forming region (age 
~ 2Myr) using the multi-object fibre spectrograph FLAMES at 
the ESO/VLT. For each target, we obtained 12 epochs, covering 
time-scales from ~ 1 week to 15 months. Rather than focusing 
on individually selected targets, our main goal was to constrain 
long-term accretion variability for typical YSOs. We use informa- 
tion from three accretion-related emission lines (Ha 6563 A, He-I 
6678 A, Ca-II 8662 A) to investigate possible accretion within these 
systems. 

The paper is structured as follows: After this introduction, we 
describe the observations and the data in Sect. [2] we discuss the 
accretion indicators and the measurements of emission lines in Sect. 
|3]and Sect.|4] Methods used to derived accretion rates are presented 
in Sect. [5] The results are discussed in Sect. [6] and conclusions are 
drawn in Sect. [7] In the Appendix we discuss the possible presence 
of veiling in the spectra of the accretors. 



2 OBSERVATIONS AND DATA PREPARATION 

The observations were carried out using the multi-object fibre 
spectrograph FLAMES/GIRAFFE at the ESO/VLT for ESO pro- 
grammes 082.C-0005(A) and 084.C-0094(A). Our targets are con- 
firmed members of the ^ 2Myr old Cha-I star forming region, se- 
lected from the census bv lLuhmanld2007l) which contains 215 ob- 
jects. Whenever possible we used the coordinates from the U CAC2 
catalogue for the fibre positioning (Zacharias et al. 2004). The 
faintest objects (spectral type later than ~M4) are not registered 
by UCAC; here we relied on the less accurate coordinates from 
2MASS. We observed one field (25 arcmin in diameter) focused on 
the area with the highest density of YSOs in the molecular cloud in 
the southern part of Cha-I with central coordinates a ll h 08 m 20.0 s , 
<5-77°36'l0" (12000). When configuring the fibres priority was 
given to objects with a) evidence for accretion from emission lines 
and b) evidence for a circumstellar disc from mid-infrared excess 
(e.g. Class II objects or classical T Tauri stars (CTTs), see Sect. 
|4j- Apart from that, no other selection criteria were used. The fi- 
nal sample of targets contained 50 objects, however twenty five 
of these objects were dropped from the analysis because no mea- 
surable stellar continuum was observed above the sky background 
emission. Two objects which have no continuum were kept as they 
both showed strong Ha emission and other spectral lines. The spec- 
tral type distribution for the objects with and without discs in this 
sample is shown in Fig.Q] 

We observed this sample 6 times in each of the two campaigns. 
The observations cover the time span from lanuary 2009 to March 
2009 and from December 2009 to April 2010. The minimum sep- 
aration between consecutive observations is ~ 1 week. The exact 
dates of the observations are given in TableQ] 

At each epoch two spectra of each source were taken with two 
high-resolution GIRAFFE grisms, HR15N and HR21, with a wave- 
length coverage of 6470 - 6790 A and 8480 - 9000 A respectively. 
These settings were selected because they cover the accretion- 
related emission lines Ha, He I, and the Ca IRT (see Sect. [3}- The 
nominal resolution for these grisms is in the range of R ~ 16000 
- 17000. For each epoch and setting we obtained a long exposure 
(800 s for HR21 and 1000 s for HR15N) and in addition a short ex- 
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Figure 1. Histogram of spectral type distribution within sample. The dashed 
line indicates all objects and the solid line indicates the objects with discs. 
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Figure 2. A spatially averaged spectrum of the sky background Ha emis- 
sion, and a single epoch spectrum from a non-accretor in the sample, 
CHXR74, which is offset by 1 for clarity. 



Table 1. Observations 


Obs. Period 1 


Time 


Obs. Period 2 


Time 


Grism 


Obs. Date 


UTC 


Obs. Date 


UTC 




04-Jan-2009 


05:10 


09-Dec-2009 


07:06 


HR21 


04-Jan-2009 


05:26 


09-Dec-2009 


07:21 


HR15 


16-Jan-2009 


04:35 


03-Jan-2010 


06:30 


HR21 


16-Jan-2009 


04:51 


03-Jan-2010 


06:46 


HR15 


03-Feb-2009 


04:40 


15-Jan-2010 


07:26 


HR21 


03-Feb-2009 


04:56 


15-Jan-2010 


07:42 


HR15 


ll-Feb-2009 


05:30 


27-Jan-2010 


04:11 


HR21 


ll-Feb-2009 


05:46 


27-Jan-2010 


04:26 


HR15 


19-Feb-2009 


05:38 


10-Feb-2010 


06:30 


HR21 


19-Feb-2009 


05:54 


10-Feb-2010 


06:45 


HR15 


Ol-Mar-2009 


02:42 


07-Mar-2010 


06:20 


HR21 


Ol-Mar-2009 


02:59 


07-Mar-2010 


06:36 


HR15 



posure (30 s) in case of saturation. As saturation was not an issue 
even for the brightest objects; we use the short exposures only for 
consistency checks. 

For our analysis, we rely on the fully processed science frames 
as provided by the ESO/GIRAFFE pipeline ('SRBS' files). The 
pipeline reduction includes bias correction, extraction, flatfielding, 
transmission correction and wavelength calibration. The spectral re- 
sponse is not corrected, because no standard stars were observed, 
i.e. absolute flux measurements are not feasible with our data set. 
Note that in the standard pipeline reduction, sky subtraction and 
scattered light subtraction does not take place. 

Spectra taken in molecular cloud regions are affected 
by emission line flux from the nebular background (e.g. 
ISicilia-Aguilar et all |2005). For fibre spectrographs this can be 
problematic, because a measurement of the background close to the 
target is usually not feasible. To be able to assess the background 
contamination, we placed a number of additional fibres on positions 
without any source (35 and 58 sky fibres for the first and second 
observation periods). Fig. [2] shows an averaged spectra from eight 
sky fibres, with a single epoch spectra from an object in the sample, 
CHXR74. Here we see that the amount the sky emission contributes 
to the Ha flux is small (here CHXR74 has a Ha EW ~ 8.6 A). 

We used two different methods to remove the sky background 
continuum levels from the object spectra. For the majority of ob- 
jects, a spatially averaged sky continuum value was derived for each 



epoch, and was removed from their spectra. To increase our spatial 
resolution and the reliability of our background subtraction, in ad- 
dition to the sky fibres, we used 18 target spectra that showed no 
continuum emission above the sky background. Using the fluctu- 
ations in our emission line measuements to test how reliable this 
sky subtraction is, we found that this method was insufficient for 
four object's shorter wavelength spectra (6470 - 6790 A), which are 
discussed in Sect[3] 



3 EMISSION LINES 

The dominant emission line in our spectra is Ha (6562.81 A) and is 
present in almost all objects. In addition, both the He I (6678.2 A) 
and Call triplet emission lines (8498, 8542, 8662.1 A) are covered 
and are present in the spectra of a few objects. 



3.1 Ha emission 

In Fig. [3] and [4] we show average line profiles for a selection of 
objects in our sample. Based on the shape of the Ha profiles, the 
objects clearly fall in two groups. Ten of them show strong, broad 
Ha emission (see Fig. |3), often asymmetric and accompanied by 
absorption features. In this case 7 out of 10 show obvious absorption 
superimposed on the emission line, 5 have blueshifted and 1 has 
redshifted absorption. Most of the remaining 15 objects have weak, 
symmetric Ha emission as seen in the first two upper panels of Fig. 
[4] A small subgroup of these objects, show broad wings in their Ha 
profiles as seen in the last two upper panels of of Fig. [4] 

Under each average profile in Fig.[3]and Fig.|4]the normalised 
variance profile is plotted, which can be used to assess the variabil- 
ity as a function of wavelength. The variance at a certain wavelength 
is given by: 



2 (A) = 



n 



(1) 



where is the flux at wavelength A for spectrum number i, I\ 
is the average flux at that wavelength and n is the total number 
of spectra for a given object jjohns &"B asri 1995). The normalised 
variance is then given by a 2 N {\) = a 2 (X)/I\. 
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The horizontal dashed line in these plots indicates the zero 
variability level, which is given by 

oir(A) = <x5v,o i^~T^j ^ 

where a % n is the normalised variance in the continuum 
dStelzer et al]|2007l ; IScholz & Javawardhanj|200^) . 

The variance profiles in Fig. [3] reveal that the variability is 
not uniform across the line profile. We can see that seven of these 
objects show multiple peaks of enhanced variations. (Two objects 
B43 and ChaHa2 are shown without variance profiles, as they have 
weak continuum emission compared to the sky background emis- 
sion). 

In all spectra we measured Ha equivalent widths (EW) and 
10% widths. The 10% width is simply the full width of the emission 
line at 10% of the peak height, and is measured in kms" . The 
equivalent width of a line is given by 

EW = E ^ir^ AA 0) 

F c 

where F c is the continuum flux and Fa is flux integrated over the 
line and the units of EW are A. We used three different integra- 
tion window sizes, depending on the emission line width (25, 15, 
or 7.5 A). The continuum was found by averaging measurements 
from two 10 A windows either side of the integration window. To 
determine the 10% width of the Ha line, we used the same contin- 
uum as for the EW to find the peak height above the continuum. In 
Table |2] we list the average and the spread in EW and 10% width 
for all objects as well as the window sizes and an estimate of the 
measurement errors. (A brief explanation of how the measurement 
errors where estimated is given in AppendixlBV 

Due to the low signal to noise ratio in some spectra, a number 
of the 10% width measurements are considered to be upper limits. 
In these cases the Ha wings are not clearly distinguishable from the 
continuum due to the noise. This does not affect the EW measure- 
ments in the same way, as to calculate the EW we integrate over the 
full line. For the objects with more measurements than upper lim- 
its, we only take the measurements in calculating the average 10% 
width given in Table [2] 

As previously mentioned in Sect|2] for four objects in the sam- 
ple, a further step in the sky subtraction was taken. When a range 
of sky continuum values from across the field of view was removed 
from these objects, it was found that their EW measurements for 
a single epoch varied by 20 - 50A, causing their derived accretion 
rate change of up to an order of magnitude. For these four objects 
(IS0143, ChaHa2, ChaHa6 and B43) we used 3 - 4 nearest neigh- 
bour fibres to estimate the local sky background continuum. Again 
by varying this local sky continuum value within the range given 
by the neighbouring sky fibres, we find the EW measurements stay 
stable for ISO 143 and ChaHa6. However objects B43 and ChaHa2 
still show very large variations in EW, and so for these objects we 
quote only 10% width measurements which are not affected by low 
continuum counts. 

In the top two panels of Fig.[5]we plot the average 10% width 
against EW for all 23 well measured objects, one panel for each 
observing period. 

Over-plotted on these two panels is the often adopted cut 
off between CTTs and Weak lined T Tauri stars (WTTs), an EW 
of lOA and a 10% width of 270 km s" 1 dWhite & Basrill2003l) . 
IBarrado v N avascues & Martini {2003) have introduced a more pre- 
cise CTTs/WTTs division based on spectral type: for K6 stars it 
is 5.9 A, for Ml stars it is 10.1 A and for M6 objects 24.1 A. The 



mean spectral type in our sample is Ml, making the 10 A a reason- 
able value to take, as the three K and two G type stars in our sample, 
either have very high EW or an EW well below the cut off for their 
spectral type. 

Eight objects stand out with EW > 20 A and 10% width 

> 350 km s -1 (two objects, B43 and ChaHa2 have 10% width 

> 270 km s" 1 but no EW measurements). They appear in the upper 
right quadrant and coincide with the objects found to have highly 
variable, mostly asymmetric line profiles (see above). Defining the 
spread as the maximum - minimum value measured during one ob- 
servation season, the spread in Ha EW and 10% width measure- 
ments are indicated by the over-plotted horizontal and vertical bars. 
Although variability is clearly present, the objects in the right upper 
quadrant have continuous strong Ha emission over both observa- 
tion periods and remain a separate group. Objects with 10% width 
< 270 kms -1 also show some variations but are not shown in the 
plot for clarity. 

The subgroup of objects mentioned above with weak Ha emis- 
sion but broad wings, are ringed in the upper plots of Fig. [5] with 
EW< 20 A and 10% width> 270 kms -1 . These objects show sig- 
nificant variability over the course of the observations; in two cases 
their 10% widths fall below 270 km s _1 in some epochs. 

It is worth noting that the changes in EW and 10% width both 
show erratic behaviour, with no recognisable pattern across the time 
period of our observations. However it is possible that the temporal 
resolution of our observations was not high enough to resolve any 
periodic or quasi-periodic variations in these quantities. 

3.2 Other emission lines 

In the spectra for which He I (6678.2 A) and Ca II (8662. 1 A) emis- 
sion lines were detected, we measured the EW of these lines. For the 
He I line we used a 6 A integration window and the average of two 
10 A windows on either side of the integration window to measure 
the continuum. For Call we use an integration window of either 
12 A or 20 A for the stronger lines, and the average of three differ- 
ent continuum measurements, two 10 A windows centred at 8605 
and 8585 A and a 20 A window centred at 8670 A. Table [3] shows 
the EW measurements for both lines. Estimates of the measurement 
errors for both lines are also given in Table[3] (See Appendix iBlfor 
an explanation of how these were derived). 

In the lower left panel of Fig. [5] we have marked the objects 
with Ca II or He I emission with squares and triangles respectively. 
It clearly shows that these emission lines are found along with 
strong and variable Ha emission. Calcium absorption increases in 
earlier types, correspondingly, the two earliest type stars in our sam- 
ple, T33A and T26, show Call in absorption. (For four epochs 
in the first observation period T33A does show weak Call emis- 
sion). Individual epoch Call EW measurements for the accretors 
are given in Table [9] Fig. [6] shows the average Call profiles for all 
10 accretors, Fig.|7]shows the average He I profiles for the five ob- 
jects that show He I emission and one which does not. 



4 ORIGIN OF EMISSION 
4.1 Accretion 

Based on the emission lines in the spectra, our sample shows a clear 
division into two groups (see Sect[3j. In this section we argue that 
10 objects in our sample are accreting from a disc for the following 
reasons: (i) Strength and width of Ha emission lines (ii) Shape of 
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Figure 3. Average and variance profiles for Ha emission feature for the 10 objects with strongest Ha emission (see Sect.[3j]. The name of the object is given 
in each frame. Vertical lines mark the central wavelength for Ha emission (6562. 8lA). The horizontal dashed lines in the variance profile plots represent the 
zero variability level allowing for signal to noise. B43 and ChaHa2 are shown without variance profiles, as their continua were not measureable. 
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Figure 4. Average Ha profiles for two objects with weak emission (left, CHXR74, ChaHa8) and for two objects with weak emission but extended wings (right, 
ESOHa560 and LM04 429). Below each is their variance profile. The horizontal dashed lines in the variance profile plots represent the zero variability level 
allowing for signal to noise. 



Ha profiles (iii) Presence of Calcium and Helium and (iv) Presence 
of circumstellar discs in the sample. 



(i) Ha 10% versus EW: Eight objects show Ha emission with 
average EW > 20 A and 10% width > 340 km s" 1 (see Fig. [5) 
(Two further objects, B43 and ChaHa2, which we classify as accre- 
tors have strong Ha emission with 10% width > 340 km s -1 , but 
do not have a continuum for measurement of the Ha EW). In the 
case of accretion we would expect the line velocities to correspond 
to the typical velocities in an accretion flow. Taking the mass range 
in our sample to be around 0.1Mq-1Mq, free fall velocities should 
be in the region of ~160 - 300 km s -1 , which would correspond to 
a range in 10% widths of ~320-600kms~ 1 . The apparent free- 
fall velocities will be smaller than this range in most cases due to 
the inclination of the stellar system with respect to us. For example 
taking an inclination of 45° this 10% width range becomes ~ 225 - 
425 km s _1 . The line widths for the 10 accretors are consistent with 
these ranges (see Fig.[5]and Table |2j. In contrast, the rotational ve- 
locities of M-type st ars in Cha-I have on average usini ~11 kms" 
dNguven et alj2009h . which cannot account for the broadening that 
we see. 

(ii) Shape of the Ha profiles: The accretors' Ha profiles are al- 
most all asymmetric, often with more than one peak and an overly- 
ing absorption feature (see Fig.[3j. Therefore the emission is likely 
to originate from a structured flow. One possible explanation of the 
asymmetry is absorption by stellar wind or m aterial in the inner disc 
jEdwards et alii 19941 ; lHartmann etai] 19941) . 

(iii) Ha Variability: The accretors show distinctively different 
variations across the emission line compared to the non-accretors. 
The variance profiles seen in Fig. [3] show multiple peaks of varia- 



tions across the emission line and are much stronger than those of 
the non-accretors (see Fig. [4j . 

(iv) Presence of Calcium and Helium emission: Strong emission 
of these l ines is thought to originate in the infall region of a magne- 
tosphere dMuzerolle et al.ll998l : |javawardhana et alj2003h . As seen 
in the lower right panel of Fig. [5] and from Table [3] 8 out of 10 
of our accretors show signs of either He I or Call emission. None 
of the non-accretors show the presence of Call or He I. Both of 
these emission lines can be produced in the chromospheres of ac- 
tive M-type objects, however it is either a t low levels compared with 
our accretors or in transitory flare eve nts jjavawardhan a et alj2006l ; 
iMontes et al.lll998l : lGizis et~ai]|2002l) . 

(v) Discs: The lower right panel of Fig. [5] shows a plot of 10% 
width versus EW of the Ha emission in our sample, where the ob- 
jects with discs are m arked with a diamon d. This is based on the 
SED classification bv lLuhman et al.l d2008l) from Spitzer photome- 
try (see Fig.[8jl. All of the objects in our sample are either Class II, 
which are considered to have inner discs, or Class III objects. All 
10 accretors are Class II objects. 

The above arguments strongly suggest that these 10 objects are 
actively accreting. This accretion is further discussed in Sect.[5]and 
[6] We define the 15 remaining objects as non-accretors (further dis- 
cussed below). There is a possibility that some are weakly accreting 
objects that fall below our detection limits. However, the distinctive 
difference in Ha emission between the accretors and non-accretors 
is strong evidence that they are two separate populations, with dif- 
ferent physical processes determining the Ha emission. 

There are four objects that do show Ha EW of ~ 10 A and 
have large Ha wings (these are ringed in FigO. Apart from these 
broad wings, they do not share the other characteristics with the ac- 
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Table 2. Data for observed sample. Spectral types from lLuhmanl (2007), SED classes from lLuhman et al] (200?). Coordinates are for epoch J2000. EW and 
10% width measurements listed are the averages from all 12 epochs. The AEW and A10% width are the maximum - minimum measurements for the 12 
epochs, unless otherwise indicated. Integration Window is the window over which the EW and 10% widths are calculated. Error estimates are given for EW 
and 10% width measurements (See Appendix [5]for more). Sky subtraction method is indicated in the second last column, M indicates that the the mean sky 
continuum value was removed for that object, NN indicates that the nearest neighbour fibre method was used to estimate the local sky continuum for removal, 
see Sect|2]for details. Note: Sky subtraction was not possible for B43 or ChaHa2 so measurements are not given. 





Spectral 


SED 




a 


5 


Ha EW 


AEW 


EiTor 


Ha 10% 


A10% 


Error 


Sky 


Integration 


Object 


Type 


Class 


fhh 


:mm:ss] 


r ' "] 


(A) 






(km s ) 






Sub. 


Window (A) 


CHXR28 


K6 


III 


11 


07:55.8 


-77 27 25.5 


0.6 


0.8 


(0.1) 


<194 




(36) 


M 


7.5 


CHXR20 


K6 


II 


11 


06:45.1 


-77 27 02.3 


0.3 * 


0.7 


(0.1) 


<255 




(24) 


M 


7.5 


T22 


M3 


III 


1 1 


06:43.5 


-77 26 34.5 


3.0 


1.3 


(0.1) 


419* 


384 


(29) 


M 


15 


IS0143 


M5 


II 


11 


08:22.4 


-77 30 27.7 


138.2* 


172.2 


(13.6) 


395 


66 


(5) 


NN 


25 


T33A 


G7 


II 


11 


08:15.2 


-77 33 53.1 


25.9 


19.7 


(0.4) 


433 * 


158 


(1) 


M 


15 


T39A 


M2 


III 


11 


09:11.6 


-77 29 12.7 


5.3 


4.2 


(0.1) 


197 


86 


(3) 


M 


15 


ChaHa2 


M5.25 


II 


11 


07:42.5 


-77 33 59.3 








349 


74 


(6) 


NN 


25 


B43 


M3.25 


II 


11 


09:47.4 


-77 26 29.0 








398 * 


68 


(5) 


NN 


25 


T45 


Ml. 25 


II 


11 


09:58.7 


-77 37 09.0 


104.4 


64.8 


(1.3) 


506 


69 


(1) 


M 


25 


ChaHa6 


M5.75 


II 


11 


08:39.5 


-77 34 16.6 


81.7 


55.7 


(3.8) 


350 


118 


(6) 


NN 


25 


ISO 126 


Ml. 25 


II 


11 


08:2.97 


-77 38 42.5 


79.9 


72.1 


(4.0) 


386 


68 


(1) 


M 


25 


ChaHa8 


M5.75 


III 


11 


07:46.1 


-77 40 08.9 


6.7 


9.1 


(1.3) 


170* 


108 


(24) 


M 


15 


ChaHo5 


M5.5 


III 


11 


08:24.1 


-77 41 47.3 


7.3* 


4.5 


(1.4) 


168* 


116 


(12) 


M 


15 


ESOHa566 


M5.75 


III 


11 


09:45.2 


-77 40 33.2 


7.5* 


3.1 


(1.0) 


t JO 


/J 


(3) 


M 


1 C 
1 J 


T30 


M2.5 


II 


11 


07:58.1 


-77 42 41.3 


59.0 


53.9 


(2.2) 


561 


213 


(1) 


M 


25 


T34 


M3.75 


III 


11 


08:16.5 


-77 44 37.2 


3.0 


2.6 


(0.2) 


194* 


78 


(11) 


M 


15 


ChaHa3 


M5.5 


III 


11 


07:52.3 


-77 36 56.9 


6.9 


10.6 


(0.5) 


222* 


146 


(12) 


M 


15 


T26 


G2 


II 


11 


07:20.7 


-77 38 07.3 


20.7 


7.6 


(0.2) 


430 


90 


(1) 


M 


25 


ESOHa560 


M4.5 


III 


11 


07:38.3 


-77 47 16.8 


10.2 


8.0 


(0.7) 


526* 


179 


(7) 


M 


15 


LM04 429 


M5.75 


III 


11 


07:24.4 


-77 43 48.9 


11.7* 


11.2 


(1.6) 


<618 




(23) 


M 


25 


CHXR22E 


M3.5 


II 


11 


07:13.3 


-77 43 49.8 


9.5* 


8.1 


(0.7) 


412* 


357 


(37) 


M 


25 


CHXR21 


M3 


III 


11 


07:11.5 


-77 46 39.4 


9.1* 


7.5 


(0.2) 


425* 


211 


(4) 


M 


25 


T31 


K8 


II 


11 


08:01.5 


-77 42 28.8 


68.9 


36.7 


(0.4) 


470 


214 


(1) 


M 


25 


CHXR76 


M4.25 


III 


11 


07:35.2 


-77 34 49.3 


12.6 


6.4 


(1.3) 


153 


58 


(2) 


M 


15 


CHXR74 


M4.25 


III 


11 


06:57.3 


-77 42 10.6 


8.6 


3.5 


(0.4) 


174* 


74 


(1) 


M 


15 


* In these cases the mean is not based on 12 measurements. See Tables[7]and[8j for 


more. 













Table 3. Emission Line Measurements for objects with either He or Ca emission in their spectra. EW measurements listed are the averages from each observation 
period, AEW is the maximum - minimum measurement in an observation period. Estimate of measurement errors are given in parentheses (See Appendix 151 
for more). 



Object 


Obs.l EW 

[A] 


Call Emission 
ACa EW Obs. 2 EW 

[A] 


AEW 


Error 


He I Emission 
Obs.l EW Obs.2EW 

[A] [A] 


Error 


ISO 143 


4.6 


5.0 


3.2 


3.5 


(0.30) 


7.7 


10.0 


(3.00) 


T33A 


2.0* 


2.5 






(0.07) 








ChaHa2 


















B43 


12.4 


15.8 


19.3 


7.0 


(0.90) 


X 






T45 


21.6 


24.7 


19.8 


16.3 


(0.05) 


1.2 


1.1 


(0.05) 


ChaHa6 












1.4 


1.3 


(0.20) 


ISO 126 


6.2 


5.6 


6.7 


4.3 


(0.10) 








T30 


11.4 


17.1 


7.4 


12.4 


(0.10) 


1.7 


1.5 


(0.20) 


T26 


















T31 


3.3 


6.2 


7.4 


8.4 


(0.01) 


0.6 


1.0 


(0.06) 



* : Ca II in emission in only four epochs for this object. 

X : Hell present in spectra, but no continuum for measurements. 



cretors. Most notably, only one object has any evidence of having 
a disc and none have Call or He I in emission. The Halpha emis- 
sion profiles are quite different to the accretors, they are centrally 
peaked and show with no overlying absorption or other structures 
within the profile. Their Ha EWs are also belo w the cut-offs for 
accretors dBarrado v Navascues & Martini 120031) . and they do not 
show the same amplitude of variations in EW (4.8 A on average) as 



we see in the accretors (36 A). This can also be seen by comparing 
the variance profiles of the non-accretors, Fig[4] with the strength 
of those of the accretors, Fig.[3] The origin of this emission remains 
unclear, however a full discussion of these objects is beyond the 
scope of this paper. 
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Figure 5. Two Upper Panels: Average Ha 10% width versus Ha EW for the first and second observation periods on the left and right respectively. The 
horizontal line is at 270 km s _1 and the vertical line is at 10 A, which are the standard limits to differentiate between WTTs and CTTs. The spread (max - min) 
in measurements for objects with the broadest Ha profiles are indicated by the over-plotted horizontal and vertical bars. The arrows pointing down indicate 
upper limits in 10% width. The ringed objects have weak but borad Ha emission and are dicussed in Sect. 14. II Lower Left Panel: Average Ha 10% width and 
EW for all observations, showing objects with Helium emission (triangles) and Calcium emis sion (squares). Lowe r Right Panel: Average Ha 10% width and 
EW for all observations where objects marked with diamonds are Class II YSOs according to Luhman et alj j200^) . 



4.2 Magnetic Activity 

The 15 non-accretors in the sample all show an average Ha 
EW of < 13 A, which is consistent with chromosphe ric emission 
from young M-type stars caused by magnetic activity dScholz et alj 
2007). They mostly show narrow, symmetric Ha profiles (see Fig. 
0, which is expected for chromospheric emission. 



5 ACCRETION RATE ESTIMATES 

The Call EW (at 8662.1 A), Ha 10% width and the Ha EW were 
used to estimate mass accretion rates. The He I emission was not 
further examined due to the fact that the line was very weak and 
noisy or else not observed at all. 

(i) Call: To deri ve M we follow the method given by 
Moh anty etai] [2005). Firstly, we determine the underlying contin- 
uum flux around the Call em i ssion line from the AMES-Dusty syn- 
thetic spectra bv lAllard et alj J200lh . which are available for a range 
of effective temperatures, gravities and metallicities. The appropri- 



ate model spect rum is selected using the effective temperature from 
Luhman (2007). In each case we used the lowest log g available, 
which for T eff < 3100K was log g = 3.5, T ef£ < 4000K was log g 
= 4 and T e a > 4000K was log g = 5. Considering a small range of 
log g values will not adversely affect our results as the continuum 
flux is only negligibly dependent on log g, and primarily determined 
by T e ff dMohantv et alj2005h . The average continuum flux (F cont ) 
was then computed from the model over the range 8600 - 8700 A. 

The Calcium flux Fcaff (per unit area on the stellar surface [ergs 
sec _1 cm -2 ]) was calculated from the EW, using the following 
equation: 

Fcaii = F cont ■ (1 + r A(CaII) ) ■ EWcatt (4) 

where rvcaii) is a measure of excess continuum emission pro- 
duced by the accretion shock relative to the photosp heric emission, 
and is given by r^(Caff) = ^f"^" a ff • Following iMohantv et alj 
(2005) we assume r A( - CaII ) ~ 0, which is a reasonable approxima- 
tion for the red part of the spectrum and at low accretion rates. If 
this is not the case, and r A (Catt) is n °t close to zero, we will have 
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Figure 6. Average and Variance profiles for Call emission in the 10 accretors. Variance profiles are only shown for objects with strong Call emission. Note 
for T33A, the average profile does not show any Call emission, but this object does show weak Call emission in four epochs. 



underestimated the accretion rates. (See Appendix [A] for a test of 
this statement). 

Using accretion rates derived from UV veiling measurements and 
detailed modeling of the H a profiles in a sample with a mass range 
of -0.02 M Q to -2 Mg^. lMohantv et alj J2005h found the follow- 
ing relation between accretion rates and Call fluxes. 



log(M) = 1.06- log (F c 



15.40 



(5) 



where the units of M are M j/r _1 . (Note|H erczeg & Hillenbrand 
(2008) found a very similar relation log(M) = 1.03 ■ log(fcaii) — 
15.2). 

(ii) Ha 10% Width: Using veiling measurements for calibration 
purposes and Ha line modeling from a few dif ferent studi e s with 
a sample mass range of 0.04 Mq to 0.8 Mq, iNatta et alj d2004l) 
found that the 10% width of the Ha emission line correlates with 
the mass accretion rate: 



log(M) = A + B ■ (Ha 10% width) 



(6) 



where A = -12.89 ± 0.03 and B = 0.0097 ± 0.0007. Units of 
M are M j/r _1 . We do not propagate the spread in the coefficients 
through our calculations, as we are interested in the variations rather 
than the absolute accretion rates. 

(iii) Ha EW: To estimate the accretion rate using the Ha EW, we 
begin by converting the Ha EW a flux per unit area on the stellar 
surface [ergs sec -1 cm~ 2 ]. As with Call and following Eq. [4] we 
use the AMES-Dusty models to compute the continuum flux in the 
range 6554 - 6572 A. Note that the AMES-Dusty models do not go 
to high enough temperatures to cover the T e ff of two objec ts T26 
and T33A. In these two cases we used the iKurucj dl993h atmo- 
sphere models to convert the Ha EW to a flux. Veiling is not taken 
into account (see Appendix lAl for further details). Using the pub- 
lished stellar luminosities and effective temperatures from[Luhrnan 
J2007I) we determine estimates for the stellar radii (R). The Ha 
fluxes are then converted to luminosities by multiplying by 4tyR 2 . 

By fitting models to accretion co ntinuum emission to objects in 
the range 0.02M Q < M < !M . lHerczeg & Hillenbrand! d2008h 
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Figure 7. Average profiles for 5 objects with He I emission and one without. 



derived a series of relations between accretion luminosity and the 
luminosity of certain emission lines associated with accretion. They 
derived the following relation for Ha: 



l0g(Lacc) = A + B • log(L Ho 



(J) 



where A = 2.0 ± 0.4 and B = 1.20 ± 0.11. L Hq , L acc are in units 
of Lq . We use this equation to determine accretion luminosities but 
we do not take into account the spread in the coefficients. Assuming 
that all the gravitational energy from the accretion is converted into 
luminosity, the accretion luminosity is re lated to the accretion rate 
as follows faercz eg & Hillenbrand 200^): 



M = 1 



Rir, 



(8) 



GM, 

Here i?» is the stellar radius, Ri n is the infall radius and M* is 



0.8 



0.6 



0.4 



0.2 



- 



r °1 5 r n 1 

r5.ai-rs.0i 

Figure 8. Colour pl ot for the sample based on Spitzer photometry from 
iLuhman et alj i2008t) . Here the accreting objects are indicated with a box. 
Four objects are not shown as not all of the 4 IRAC magnitudes were avail- 
able for them (T26, T33A, ESOHa560 and CHXR21). 



Table 4. Parameters used in the derivation of accretion rates. Fluxes are in 



units of [erf 


;s sec 1 cm 2 < 


4.- 1 ]xio fi 






Object 


Fcont(Ha) 


Fcont(Ca) 


M» 










[M ] 


[R©] 


IS0143 


1.369 


3.600 


0.175 


1.02 


T33A 


58.89* 




2.5 


1.77 


ChaHa2 






0.15 


1.17 


B43 




5.731 


0.35 


1.38 


T45 


4.383 


8.073 


0.75 


2.29 


ChaHa6 


0.879 


2.797 


0.11 


0.87 


IS0126 


4.383 


8.073 


0.75 


1.43 


T30 


3.055 


6.467 


0.50 


1.07 


T26 


70.16* 




2.5 


3.39 


T31 


10.295 


11.567 


1.15 


3.71 



Estimates from Kurucz (1993) models. 



the stellar mass. Following Herczea & Hillenbrand! (2008), we ap- 
proximate the factor ( 1 — ~ 1.25, assuming Ri n ~ 5R* 
bullbring etai]|l998l) . 

To get an estimate of the stellar mass, we compare the pub- 
lished effective temperatures of our objects ( lLuhmarj [2007) with 
a 2Myr isochrone from the BCAH98 evolutionary models by 
lBaraffeetai1dl998h . Two objects in our sample have temperatures 
that are higher than those covered by the Baraffe models (T26 
and T 33A); for these objects we e stimate M ~ 2 .5 Mp based 
on thelry Antona & Mazzitell| dl997h isochrones. The lBaraffe et alj 
( 1998) models have bee n found to systematical l y ove restimate 
masses compared to the ID' Anto na & Mazzitellij j 19971) models 
(Muzerolle et al. 2003). Therefore we have underestimated the Ha 
EW derived accretion rates of T26 and T33A compared to those of 
the other accretors. This means that although these accretion rates 
are not fully comparable, the relative spread in accretion rates are. 

The time-average accretion rates and the spread (max - min) 
in accretion rates for all three diagnostic are given in Tables [5] and 
[6] The parameters such as F cont , and R, that were used in the 
derivation of these accretion rates are given in Table [4] 
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Table 5. Derived time-averaged accretion rates for the sample in Table|4] Following each average is the max - min spread in accretion rates over the 12 epochs 
of observations. An estimate of measurement errors are give in parentheses. Units: Log(M0yr~ 1 ). A 2a upper limit is given for the object with no discernible 
Call emission in it's spectrum. No Call accretion rate estimates are given for the 3 objects in our sample with Call in absorption. For B43 and ChaHa2 no 
Ha EW estimates are given as they have no measurable continuum at that wavelength. 



Object 


Log(M) 


Error 


Log(M) 


Error 


Log(M) 


Error 




Ca EW 




Ha EW 




Ha 10% 




IS0143 


-o.y/ ± u.o/ 


(0.07) 


o qt _i_ n in 
-o.oz ± o.zy 


(0.10) 


-y.oo ± u.o4 


(0.10) 


T33A 






-8.09 ± 0.42 


(0.01) 


-8.68 ± 1.53 


(0.01) 


ChaHa2 










-9.54 ± 0.68 


(0.12) 


B43 


-8.39 ± 1.04 


(0.09) 






-9.02 ± 0.70 


(0.09) 


T45 


-7.85 ± 0.68 


(0.003) 


-7.82 ± 0.34 


(0.01) 


-7.98 ± 0.66 


(0.01) 


ChaHa6 


<-9.16f 


(0.08) 


-9.37 ± 0.40 


(0.05) 


-9.49 ±1.14 


(0.12) 


ISO 126 


-8.39 ± 0.73 


(0.04) 


-8.65 ± 0.46 


(0.06) 


-9.14 ±0.66 


(0.01) 


T30 


-8.39 ± 1.11 


(0.03) 


-9.24 ± 0.59 


(0.04) 


-7.45 ± 2.07 


(0.02) 


T26 






-7.16 ±0.19 


(0.01) 


-8.71 ±0.88 


(0.01) 


T31 


-8.20 ± 0.76 


(0.03) 


-7.06 ± 0.27 


(0.01) 


-8.32 ± 2.08 


(0.01) 


Av. Spread 


± 0.83 




± 0.37 




± 1.11 





t: 2a upper limit 



Table 6. Derived average accretion rates shown separately for observation periods 1 and 2. Units: Log(MQj/r — 1 ). The 2<r upper limits are given for the 
object with no discernible Call emission in it's spectrum. For B43 and ChaHa2 no Ha EW estimates are given as they have no measurable continuum at that 
wavelength. No Call accretion rate estimates are given for the 3 objects in our sample with Call in absorption. 





Observation Period 1 


Observation Period 2 






Log(M) 


Log(M) 


Log(M) 


Log(M) 


Log(M) 


Log(M) 




CaEW 


Ha EW 


Ha 10% 


Ca EW 


Ha EW 


Ha 10% 


ISO 143 


-8.89 ± 0.58 


-8.80 ±0.18 


-9.06 ± 0.50 


-9.04 ± 0.44 


-8.96 ± 0.01 


-9.05 ± 0.64 


T33A 




-8.02 ± 0.25 


-9.00 ± 0.52 




-8.20 ±0.34 


-8.31 ± 1.10 


ChaHa2 






-9.65 ± 0.42 






-9.45 ± 0.47 


B43 


-8.24 ± 0.64 




-8.91 ± 0.23 


-8.60 ± 0.50 




-9.26 ± 0.65 


T45 


-7.78 ± 0.68 


-7.82 ± 0.34 


-7.84 ±0.51 


-7.88 ± 0.54 


-7.83 ±0.14 


-8.11 ±0.39 


ChaHa6 


<-9.26 f 


-9.36 ± 0.40 


-9.44 ± 1.12 


<-9.11 f 


-9.38 ± 0.24 


-9.55 ± 1.05 


ISO 126 


-8.45 ± 0.64 


-8.70 ± 0.30 


-9.14 ± 0.66 


-8.35 ± 0.64 


-8.61 ±0.36 


-9.13 ±0.33 


T30 


-8.28 ± 1.10 


-9.28 ± 0.59 


-7.07 ± 0.94 


-8.49 ±1.10 


-9.23 ±0.31 


-7.83 ±2.01 


T26 




-7.13 ±0.10 


-8.81 ±0.69 




-7.20 ±0.11 


-8.62 ± 0.64 


T31 


-8.36 ± 0.50 


-7.10 ±0.16 


-8.45 ± 1.31 


-8.06 ± 0.59 


-7.02 ± 0.24 


-8.20 ± 2.08 


Av. Spread 


±0.69 


±0.29 


±0.69 


±0.54 


±0.22 


±0.94 



f : 2a upper limit 



6 DISCUSSION 

6.1 Comparing accretion rate indicators 

In Fig.[9]a graphical comparison of the three accretion rate estimates 
can be seen, with the 1-1 relation indicated by the dotted line. The 
over-plotted horizontal and vertical bars do not show the errors in 
measurements, which are substantially smaller (see Table[5] These 
bars denote the max - min spread in derived accretion rates from 
the 12 epochs of observations. The data points with arrows mark 
the 2a upper limit for the object with no discernible calcium in 
its spectrum. Three further objects had Call in absorption, and so 
we have no Call measurements for these three objects. Taking into 
account the spread in derived accretion rates of individual objects, 
there is a reasonable agreement between the three accretion rate 
estimates. Some more additional trends can be seen by comparing 
the average accretion rates and the amplitude of variations. 

On average the Ca EW and the Ha EW give the lowest vari- 
ations in accretion rates over the course of our observations (see 
top panel of Fig. [9j. For a given object the Ca EW estimates show 
greater variations than the Ha EW estimates do. This could possi- 



bly be explained by the presence of telluric emission lines in the 
wavelength domain of Call, which could contaminate the EWs 
dOsterbrock et alj[l996t) . The telluric lines are drowned out by the 
Call emission when it is strong (e.g. in the case of T31 and T45, 
see Fig. [6J, but they can clearly be seen in objects with weak Call 
emission. These telluric lines will change over the course of our ob- 
servations, and increase the variability that we see in the Ca EW. 
This does not account for all of the variability though, as some of 
our objects with the largest spread in Ca EW (e.g. T45 and T31) 
show no sign of telluric lines in their spectra. 

The variations in individual object's accretion rates are larger 
when derived from the Ha 10% width (~1.11 dex) than they are 
when derived from the other two signatures. (Note: 1 dex means 
there is an order of magnitude spread in measurements). These re- 
sults suggest that the Ha 10% width does not give reliable estimates 
of average accretion rates, especially with single epoch observa- 
tions. 

The fact that the amplitude of variations for individual accre- 
tion rates derived from Ha 10% width is on average greater than 
what it is for the other two accretion rate estimates implies there 
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Figure 9. (a) Calculated accretion rates from Ca EW and Ha EW. Over- 
plotted horizontal and vertical bars indicate spread in measurements for all 
1 2 epochs. The object with a downwards pointing arrow is the 2a upper limit 
for the Ca EW accretion rates as this object has no discernible calcium in its 
spectrum, (b) Calculated accretion rates from Ha EW and Ha 10% width, 
(c) Calculated accretion rates from Ca EW and Ha 10% width, (see Sect. 
[5J. The dotted line in these three plots indicates the 1-1 relation between the 
two indicators. 



is another contributor to the 10% width changes. It has been shown 
that ou tflows do contribute to the Ha emission in the wings in some 
cases JWhelan et ai]l2009l : lAlencar et al.ll200lh . Looking at the av- 
erage and variance profiles in Fig [5] we can see that much of the 
variability in the Ha line is found in the wings. For the calcium 
emission lines we see a different behaviour: we found that the Ca 



line emission varies all across the profile and not just in the wings 
(see Fig.|6j. 

Most of the accretors show blue- or redshifted absorption in 
their Ha profiles. The blue absorption is unlikely to be associated 
with accretion and more likely to be a consequence of a stellar wind. 
Over the course of our observations we see these absorption features 
change in strength and wavelength. This will contribute to the vari- 
ability in the 10% width, but it could also reduce the peak height 
of the Ha emission line, which would lead us to overestimate the 
10% width and the accretion rate. Hence, the absorption features 
could explain the fact that accretion rates derived from the 10% 
width are more variable. Furthermore, in a number of profiles there 
is evidence of high velocity absorption in the wings, which the 10% 
width measurements are far more sensitive to than the EW mea- 
surements. From these results we consider the Ca EW and Ha EW 
estimates to reflect the accretion rates more accurately. In the fol- 
lowing we focus the discussion on these two quantities and do not 
consider the Ha 10% width estimates. 

The variability cannot entirely be due to our measurement er- 
rors. Estimate of the measurements errors are given in Tables[2]and 
[3] These have been carried through to the accretion rate estimates 
and are given in Table [5] In all cases these errors are less than the 
variations in accretion rates. 



6.2 Spread in M oc A/* relation 

The average spread in accretion rates is ~ 0.83 dex for the values 
derived from the Ca EW and ~ 0.37 dex for the Ha EW. We can 
take this as the maximum spread in accretion rates for young ob- 
jects (~ 2 Myr) over a time period of 15 months. This spread is 
an upper limit as it is likely that emission from the wind and the 
chromosphere both contribute to these lines. 

This spread is relevant in the context of the M — M, relation 
(see Sect. QJ. As the size of our sample is small, we do not aim to 
derive such a relation from our dataset. However, it is found that 
when these two stellar properties are plotted against each other the 
dominant feature is the large spread in accretion rates for any given 
mass. Our multi-epoch observations are designed to probe this sig- 
nifica nt attribute of t he M — M, relation. 

iMohantv et alj {2005) derived a correlation M oc M* for a 
sample with a mass range of 0.15 Mq to ~ 2 Mq and age range 
from > 1 Myr to 10 Myr. Apart from the correlation, there is a ± 1 .5 
order of magnitu de spread in accretion rates at any given mass. 
iNatta et ail J2006h found a similar trend for a sample of (0.03 - 
3Mq ) in Ophiuchus, which has an age of about (>0.5 Myr - 1 Myr). 
For this single region there is still a spread of ~2 dex at least for 
any given mass. Our results show that this spread cannot be due to 
variations of individual objects on the time-scale of about a year. 
This statement was also checked using the standard deviation of the 
accretion rates within the sample, which on average for the Ha EW 
derived accretion rates is 0.16, for the Ha 10% width derived ac- 
cretion rates it is 0.61 and for the Call EW derived accretion rate it 
is 0.25. 

To illustrate this, in Fig.[l0]we have plotted the mean Ha EW 
and Ca EW derived accretion rates versus mass for our sample. The 
over-plotted vertical bars in this plot indicate the spread in mea- 
surements over the 12 epochs of observations. The dotted line on 
both o f these plots is the M-M* relation found by Moha nty et alj 
J2005h . Also plotted as a dashed line is the linear fit derived from 
our data set using standard linear regression. The relation derived 
from our sample is much shallower for the Ca EW derived accre- 
tion rates than the Ha EW derived accretion rate. Significantly, our 
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plot shows that the spread around the relation is much larger than 
any spread we see in accretion rates over the observations. 

Several alternative explanations for the scatter in M have been 
suggested in the literature. It could be accounted for by evolutionary 
differences between sou r ces, a s M is expected to change with time. 
However as Natta et alj J2004h still find large spreads in accretion 
rates within single star formation regions, it suggests that evolu- 
tionary effects are unlikely to explain the full spread. Differences 
in initial conditions or environmental conditions between objects 
could also result in a variety of accretion rates iDullemond et alj 
20061). ?or instance, as a resu lt of different initial disc conditions 
Alexander & Armitagd d2006h have predicted a 2dex spread in 
accretion rates for stars in the mass range 0.02-0.2 Mq and a 
2. 8 dex range for masse s 0.4- 4.0 Mq. Observations such as those 
by Rigliaco et alj d201lh have provided evidence for a mass depen- 
dent evolution in the accretion rates, with very low mass stars and 
brown dwarfs evolving faster, steepening the slope of the M oc M* 
relation at the low mass end. Furthermore, we cannot rule out 
that a significant portion of the scatter in the M - M* plots could 
be due to uncertainties in the determination of the stellar masses 
dAntoniucci et alj201 lh . stellar radius, infall radius or from the veil- 
ing calibrations themselves. 

The exponent in this relation is also poorly defined and 
there has been a range of v alues found, fr om M-Ml' 1 
tearentsen et al.l201 lh. M-Ml 8 dNatta et al J2"006h and up to M - 
M'i A dFang et al. 2009). There are also suggestions that a close to 
qua dratic relation is a resul t of limited sample siz e and detectabil- 
ity dClarke & Pringlejl2006l : lBarentsen et al]|201ll) . although some 
relation between mass and accretion rate is expected. However all 
these relations have a large spread of at least 1 order of magnitude 
in accretion rates for a given mass. 



6.3 Time-scales of Variations 

In the first observation period, the mean amplitude of variations 
in M is ~ 0.69 dex when derived from Call and ~ 0.29 when 
derived from Ha EW. For the second observation period they are 
~ 0.54 and ~ 0.22 dex respectively. The total variations over 15 
months (0.83 and 0.37 dex, see above) are not significantly differ- 
ent. The slight increase in the amplitude of variations in the full 
sample shows there is some variability on the longer time-scales. 
However, the fact that it does not increase by much means the dom- 
inant variations occur in a single observation period i.e. it is shorter 
than 8 weeks. 

Fig. QT] shows the average difference in accretion rates on dif- 
ferent time-scales within our sample. The accretion rates are those 
derived from the Ha EW (B43 and ChaHa2 are not included due 
to the lack of continuum around Ha) . In order to cover all the 
possible time-scales in our sample, every accretion rate estimate is 
compared to every other accretion rate estimate for that object. The 
differences in accretion rates (M1-M2) are then placed into three 
different time bins, depending on the time difference between mea- 
surements, (i) T < 25 d(ays) (ii) 40 d < T < 90 d (iii) T > 1 yr. The 
mean accretion rate difference in each bin is then plotted for each 
object. For two objects the largest variations occur on the shortest 
timescales. Five more objects show a roughly even distribution over 
the three time bins, and one object shows the largest variations on 
the longest timescales, T33A. This could be a true increase in accre- 
tion rate, or it could also come from the fact that we have incomplete 
sampling of the variations. Although long term variations are seen 
in the sample, for all, excluding T33A, the amplitude of variations 
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Figure 10. (a) Ha EW derived accretion rates versus object mass, (b) Ca 
EW derived accretion rates versus object mass. T he dotted line represents 
the M-Ml relation found bv lMohantv et alJfeOOsl) . The dashed line repre- 
sents the lineal' fit derived from our measurements. The over-plotted vertical 
bars indicate max - min spread in measurements. Upper limits are indicated 
by arrows. 



have reached a maximum or are within 70% of the maximum after 

25 d ays. 

iNguven et al.1 d2009l) observed spectroscopically a number of 
low-mass pre-main sequence stars over a series of time-scales of 
days, weeks and months, with the majority of the observations tak- 
ing place on shorter time-scales than our observations. They found 
that the amplitude of variations increases on time-scales from hours 
to several days, after which it saturates. Our observations are able to 
test if there are any larger variations when the period of observations 
is extended. Our results suggest that observations on time-scales of 
a few of weeks are sufficient to limit the majority of the extent of 
accretion rate variations in typical young stars over the time period 
of ~ 1 year. 

A previous long term photometric monitoring program also 
found low levels of variability in a large sample 72 CTTs over the 
course of ~20yrs dGrankin et al.ll2007f) . Only a fraction of their 
sample showed significant changes over time-scale of months to 
years, suggesting that T Tauri photometric variability is dominated 
by short term variations. 

A major proportion of the variations occur on the time-scale 
of 8-25 days, which is on the sam e order as the rota tion period 
of young stellar objects, 1-10 days dHerbst et alj|2007t) . Note that 
the rotation periods in our sample are unknown. This suggests 
the origin of the variations is close to the surface of the star, 
possibly within the co-rotation radius. If rotation is the underly- 
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ing mechanism that causes the observed accretion variability, we 
could be tracing structural asymmetries in the accretion flow rather 
than actual accretion rate changes. These time-sc ales could also 
be in agreement with mag netic reconnection events l lGoodson et alJ 
ll998l ; lHavashi et alll996h . but we do not have the temporal resolu- 
tion to distinguish what is occurring on these shorter time-scales. 

As can be seen in Fig. [5] our accreting objects are steadily 
accreting through all epochs in our observations. We do not see 
a radical change in the emission lines of the non-accretors over 
the course of our observations, they always have much weaker Ha 
emission than the accretors. This shows that if they are accreting it 
is below our detection levels, and never accrete at the same levels 
as the 10 accretors in the sample. Extreme variations in accretion 
or episodic accretion is thou ght to be necessary for YSOs to gain 
the masses that we observe fevansetalj 120091) . There have now 
been a number of observations of FU Ori events where an object 
will suddenly jump from having an accretion rate of 10 -7 MQyr -1 
to 10 _4 MQyr _1 in a very short amo unt of time, and then deca y 
slowly over the course of 20-100 yrs ( Hart mann & Kenvod[l996l) . 
However these are rare events, only a handful have been found, and 
they are expected to occur at the early stages of stellar evolution 
when the infall is still occurring to the disc. Slightly less dramatic 
accretion changes occur in EXors type objects, which can s how 1 - 
3 mag nitude flare ups every few years for a period of weeks (Herbia 
120071) . 

Highly variable, sporadic accretion rate changes have also 
bee n found in older obi ects such as in the 8 Myr o l d rj Cha clus- 
ter jMurphvetal.ll201ll : Ijavawardhana et alj[20o3) . iMurphv et all 
suggest that such sporadic accretion could occur at 'the crit- 
ical age range of 5-10 Myr when inner discs are being cleared 
and giant planet formation takes place'. The T Tauri star, T 
Chamaeleontis, shows very large cha nges in Ha EW, incr easing 
from A to 30 A and falling again, Schi sano et alJ (2009) infer 
these variations to be a result of an evolved, clumpy disc. Further- 
more Ijayw^dhanaetal] {2006) have also found evidence to sug- 
gest that inner disc clearing and/or grain growth might indeed affect 
two accretors in rj Cha. 

Our sample does not show signs of episodic or sporadic ac- 
cretion like the cases discussed above. The fact that at 2 Myr the 
younger stars in Cha I are co ntinuously accreting, supports the re- 
sults of IMurphv et alJ d201ll) that the sporadic accretion is likely 
related to the older age of the population and the evolutionary stage 
of the disc. 



7 CONCLUSIONS 

We present in this paper a study designed to investigate the accre- 
tion variability of YSOs based on multi-epoch optical spectra of 
25 targets in the Chameleon I region. Our observations took place 
over a total of 15 months. They were split into two separate peri- 
ods, of 8 weeks each. In each of these periods 6 observations were 
made, evenly spread across the 8 weeks. Our high resolution spec- 
tra from the FLAMES/GIRAFFE spectrograph on the ESO-VLT 
covered the Ha (6562.81 A), Hel (6678.2 A) and Call (8662.1 A) 
emission lines which we used to identify 10 accretors in our sam- 
ple. We derived the accretion rates of these objects from the Ha 
EW, Ha 10% width and the Ca EW. 

The accretion rate variability in our sample reached a maxi- 
mum or was within 70% of the maximum within ~ 8 - 25 days, for 
7 out of 8 accretors (based on Ha EW derived accretion rate, and 
hence excluding two accretors B43 and ChaHa2). This result com- 
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Figure 11. Average spread in accretion rate over different time-scales for ac- 
cretors in the sample, excluding B43 and ChaHa2, as these accretion rates 
were derived from Ha EW. Three different time bins were used (i)8 d(ays) 
< T < 25 d (ii) 40 d < T < 90 d (iii) T > 1 yr. Seven out of the eight ac- 
cretors here show very little increase in variations across the three time bins, 
demonstrating that the variations are dominated by short term variability. 



pleme nts those from previo us studies which covered shorter time- 
scales dNguven et aD l2009). This period of variations roughly cor- 
responds to the typical rotation period of these types of objects, 
suggesting the rotation could be modulating the observed accre- 
tion in these systems. Variability on these time-scales could be ex- 
plained by several scenarios described in Sect. Q] One such model 
is the asymmetric accretion flow hypothesis due to misaligned ro- 
tation and field axes or complex field topologies. In this case we 
would be probing the structural asymmetries in the accretion flow 
as different parts of the flow come into view. Differential rotation 
between the disc and the star has been shown to produce magnetic 
field line expansion and reconnection over the course of a few days 
dHavashi et al.ll996l) . and could also account for the variations. Fur- 
ther studies with greater temporal resolution are needed to discrim- 
inate between such hypotheses. 

When we considered the amplitude of variations in accretion 
rates from all three accretion rate estimators, we found that Ha 
10% width accretion rate estimates of individual objects show much 
greater variability than the other two estimators. For this reasons we 
do not consider the Ha 10% width as a good quantitative accretion 
rate estimator. 

Considering the full observation period we found the average 
spread in accretion rate to be 0.83 dex for Ca EW estimates and 
0.37 dex for the Ha EW estimates. These variations demonstrate 
that ove r < 15 months, the ac cretion rate changes are 0.37 - 0.83 dex 
or less. I Mohanty et alJ d2005l) . using the Call (8662A) emission 
line, found the spread for a single mass in the M-Af* relation to 
be 1-2.5 orders of magnitude. Bearing in mind that our Call sam- 
ple is small, our results show that the spread in the M-M* relation 
is probably not due to the variability of individual objects on the 
time-scales of years. They also suggest that this spread is not en- 
tirely due to evolutionary differences within a sample, but is more 
likely due to environment or initial conditions or uncertainties in 
the derived parameters. 
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Table 7. List of Ha EW values for all observations. Here A EW is the Max EW - Min EW measurement. Also given are the Mean EW and RMS values which are calculated over the entire sample. B43 and ChaHo2 
Ha EW measurements are not included here as they have no measurable continuum at that wavelength. 



Object 


EW1 


EW2 


EW3 


EW4 


EW5 


EW6 


EW7 


EW8 


EW9 


EW10 


EW11 


EW12 


Mean 


Max 


Min 


AEW 


RMS 


CHXR28 


0.4 


0.7 


0.6 


0.5 


0.7 


0.4 


0.6 


0.6 


1.2 


0.4 


0.5 


0.6 


0.6 


1.2 


0.4 


0.8 


0.2 


CHXR20 


0.1 


- 


0.7 


0.6 


0.2 


0.7 


0.8 


- 


0.2 


0.2 


- 


0.5 


0.3 


0.7 


0.1 


0.7 


0.3 


T22 


3.1 


2.2 


3.1 


3.4 


3.3 


3.2 


3.2 


2.7 


- 


3.1 


2.7 


3.5 


3.0 


3.5 


2.2 


1.3 


0.4 


ISO 143 


172.2 


129.2 


122.5 


163.4 


- 


160.6 


119.6 


- 


- 


- 


- 


9.6 


138.2 


172.2 


99.7 


72.6 


27.2 


T33A 


35.7 


22.1 


26.0 


30.9 


30.8 


34.8 


17.0 


30.6 


21.1 


19.2 


3.1 


6.1 


25.9 


35.7 


16.1 


19.7 


9.4 


T39A 


4.5 


4.8 


5.1 


8.6 


4.7 


4.4 


7.0 


4.7 


4.8 


4.7 


5.3 


4.9 


5.3 


8.6 


4.4 


4.2 


1.3 


T45 


94.2 


70.3 


93.7 


121.7 


116.6 


135.1 


107.6 


97.9 


101.1 


124.8 


94.6 


5.1 


104.4 


135.1 


70.3 


64.8 


17.7 


ChaHa6 


97.5 


83.2 


74.4 


103.4 


47.7 


93.9 


78.2 


94.4 


63.3 


100.1 


- 


3.1 


81.7 


103.4 


47.7 


55.7 


18.1 


ISO 126 


50.2 


81.4 


82.9 


89.8 


60.7 


70.7 


76.2 


79.2 


122.4 


61.7 


104.8 


8.8 


79.9 


122.4 


50.2 


72.1 


19.6 


ChaHa8 


6.4 


9.0 


- 


- 


2.7 


4.6 


- 


- 


- 


- 






5.6 


9.0 


2.7 


6.3 


2.7 


ChaHa5 


10.2 


- 


8.2 


7.1 


6.7 


5.8 


- 


- 


- 


5.8 


- 


- 


7.3 


10.2 


5.8 


4.5 


1.7 


ESOHa566 


9.3 


- 


- 


7.1 


6.2 


7.5 


- 


- 


- 


- 


- 


- 


7.5 


9.3 


6.2 


3.1 


1.3 


T30 


79.5 


52.8 


43.8 


56.0 


79.8 


26.0 


53.8 


41.3 


- 


74.7 


75.3 


6.2 


59.0 


79.8 


26.0 


53.9 


17.7 


T34 


3.0 


3.8 


3.1 


2.0 


1.8 


2.6 


3.0 


3.1 




2.6 


4.4 


4.0 


3.0 


4.4 


1.8 


2.6 


0.8 


ChaHa3 


5.4 


10.5 


6.4 


4.4 


2.1 


5.5 


5.7 


4.0 


12.2 


2.7 


12.7 


0.8 


6.9 


12.7 


2.1 


10.6 


3.7 


T26 


20.4 


24.7 


22.0 


20.3 


22.3 


22.6 


17.2 


17.1 


20.9 


19.4 


20.2 


1.2 


20.7 


24.7 


17.1 


7.6 


2.2 


ESOHa560 


11.4 


14.6 


13.2 


9.9 


7.9 


8.9 


6.6 


8.3 


11.4 


9.0 


11.4 


0.1 


10.2 


14.6 


6.6 


8.0 


2.3 


LM04 429 


13.9 


17.8 


14.4 


13.4 


6.6 


9.7 


12.5 






6.8 


11.8 


0.3 


11.7 


17.8 


6.6 


11.2 


3.5 


CHXR22E 


11.8 


9.4 


11.1 


8.7 


4.6 


9.1 


5.8 


11.0 




9.2 


12.7 


1.4 


9.5 


12.7 


4.6 


8.1 


2.5 


CHXR21 


12.3 


7.3 


12.2 


10.0 


4.8 


9.2 


5.7 


11.9 




10.2 


11.2 


5.6 


9.1 


12.3 


4.8 


7.5 


2.8 


T31 


74.7 


59.1 


62.2 


55.1 


67.4 


62.1 


79.2 


72.9 


91.7 


58.0 


60.6 


3.4 


68.9 


91.7 


55.1 


36.7 


11.5 


CHXR76 


12.2 


12.7 


16.7 


11.7 


15.0 


14.6 


11.2 


10.8 


10.8 


10.6 


14.3 


0.3 


12.6 


16.7 


10.3 


6.4 


2.1 


CHXR74 


7.4 




10.1 


9.3 


7.3 


8.6 


10.0 


8.2 


9.9 


8.5 


8.5 


6.6 


8.6 


10.1 


6.6 


3.5 


1.2 



5 
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Table 8. List of Ha 10% width values for all observations. Here A10% is the Max - Min 10% width measurement. Also given are the Mean 10% width and RMS values which are calculated over the entire sample. OS 



Object 



10% width [km s _1 ] 





1 


2 


3 


4 


5 


6 


7 


8 


9 


10 


1 1 


12 


Mean 


Max 


Min 


A 10% 


RMS 


CHXR28 


<115 


<219 


<130 


<121 


<258 


<122 


336 


257 


<252 


<206 


<266 


<252 


<194 


- 


- 


- 


- 


CHXR20 


<119 


<324 


<216 


<175 


<336 


<175 


<338 


<293 


<221 


<219 


<298 


<339 


<255 


- 


- 


- 


- 


T22 


393 


<492 


<370 


<362 


384 


<400 


407 


438 


<145 


415 


416 


414 


418* 


439 


384 


54 


17 


IS0143 


381 


369 


391 


387 


421 


416 


409 


395 


368 


435 


390 


374 


395 


435 


368 


66 


21 


T33A 


430 


391 


392 


432 


379 


381 


430 


537 


475 


424 


<380 


492 


433* 


537 


379 


158 


50 


T39A 


174 


212 


177 


261 


182 


187 


213 


188 


193 


188 


178 


204 


196 


261 


174 


86 


24 


T45 


526 


486 


521 


527 


539 


519 


496 


510 


477 


502 


470 


499 


506 


539 


470 


68 


21 


ChaHa6 


361 


300 


360 


381 


313 


417 


308 


385 


308 


407 


358 


299 


350 


417 


299 


117 


42 


ISO 126 


427 


403 


380 


358 


379 


370 


412 


379 


378 


385 


380 


390 


386 


427 


358 


68 


18 


ChaHa8 


203 


178 


195 


192 


<221 


202 


159 


186 


95 


187 


161 


112 


170* 


203 


95 


108 


36 


ChaHa5 


154 


131 


195 


155 


213 


154 


149 


154 


163 


133 


248 


<115 


168* 


248 


131 


116 


36 


ESOHa566 


164 


<133 


174 


146 


124 


154 


144 


198 


- 


162 


176 


140 


158* 


198 


124 


74 


21 


T30 


612 


623 


643 


610 


546 


564 


429 


509 


522 


637 


534 


499 


561 


643 


429 


213 


65 


T34 


176 


<254 


188 


214 


195 


200 


149 


201 


154 


226 


223 


209 


194* 


226 


149 


77 


26 


ChaHa3 


265 


217 


267 


269 


<245 


122 


218 


241 


166 


207 


247 


216 


222* 


269 


122 


146 


45 


T26 


408 


452 


437 


439 


381 


405 


405 


451 


412 


471 


435 


465 


430 


471 


381 


90 


27 


nISUrtajoU 


517 


645 


578 


49 / 


516 


err 
555 


405 


478 


<535 


5 1 5 


49 1 


533 


5ZO 


645 


405 


179 


5 1 


LM04 429 


<721 


701 


<635 


<633 


<463 


<726 


<786 


<731 


<465 


<574 


<564 


<497 


<618 










CHXR22E 


572 


393 


539 


427 


214 


469 


274 


542 


<168 


410 


444 


270 


412* 


572 


214 


357 


199 


CHXR21 


494 


371 


485 


420 


378 


436 


284 


492 


<327 


479 


443 


393 


425* 


494 


211 


494 


283 


T31 


478 


436 


477 


443 


523 


388 


514 


528 


495 


487 


330 


544 


470 


544 


330 


214 


62 


CHXR76 


151 


148 


149 


171 


167 


156 


146 


148 


125 


148 


142 


184 


153 


184 


125 


58 


15 


CHXR74 


186 




155 


185 


187 


207 


162 


133 


201 


150 


177 


170 


174 


207 


133 


74 


22 


B43 


403 




422 


399 




416 




417 




386 


350 




399 


422 


350 


72 


25 


ChaHa2 


317 


339 


336 


361 


326 


324 


388 


346 


339 


356 


342 


357 


344 


388 


317 


70 


19 
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£ 



TO 

TO 
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>3 
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I 

<3 



Indicates the cases where the mean 10% width is not based on all 12 epochs due to some measurements being upper limits. 
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Table 9. All Call EW measurements for the accretors with discernible Call emission. AEW is the Max - Min measurements. The RMS and the mean are calculated across the entire observation period. 



Object 


EW1 


EW2 


EW3 


EW4 


EW5 


EW6 


EW7 


EW8 


EW9 


EW10 


EW11 


EW12 


Mean 


Max 


Min 


AEW 


RMS 


IS0143 


4.0 


2.6 


2.0 


4.2 


4.8 


7.1 


3.6 


1.6 


4.6 


4.1 


2.0 


2.3 


3.6 


7.1 


1.6 


5.5 


1.6 


T33A+B 


2.3 


1.7 


1.1 




- 


1.4 


- 








- 




1.6 


2.3 


1.1 


1.2 


1.0 


ChaHa2 




































B43 


14.4 




21.5 


6.0 




5.4 




6.5 




6.1 


2.3 




8.9 


21.5 


2.3 


19.2 


6.5 


T45 


16.3 


7.3 


12.1 


22.3 


31.9 


26.9 


21.7 


21.6 


23.9 


17.2 


10.2 


7.3 


18.2 


32.0 


7.3 


24.6 


7.8 


ChaHa6 




































ISO 126 


2.7 


6.5 


5.1 


6.4 


1.8 


7.3 


4.9 


4.3 


9.9 


4.9 


7.0 


4.6 


5.5 


9.9 


1.8 


8.0 


2.3 


T30 


18.7 


5.8 


8.2 


11.3 


15.1 


1.7 


3.7 


2.2 


19.3 


15.4 


5.3 


3.0 


9.2 


19.3 


1.7 


17.6 


6.8 


T26 




































T31 


6.9 


2.7 


3.1 


2.4 


5.7 




12.4 


10.8 


6.0 


3.5 


5.2 


9.8 


6.2 


12.4 


2.4 


10.1 


3.0 
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-0.4 



-0.6 



~0.f 




Table Al. Accretion rate estimates derived from Ha EW with and without 
veiling correction for four accretors, and calculated average veiling factor, 
r. Units of M are [Mg^r -1 ] 



IS0126 



10 



Epoch Number 



Figure Al. Li EW measurements of three accretors and a non-accretor 
(T39A) used as a reference. Note: The measurements for the fourth accre- 
tor and its reference non-accretor are not shown in this plot, as they are of 
different spectral types. 



APPENDIX A: TEST FOR VEILING 

As mentioned in Sect. [T] accretion processes also emit continuum 
emission. It fills in and 'veils' any photospheric absorption lines, 
making these lines appear shallower. One method of estimating the 
accretion rate is to compare a veiled spectrum of an accreting star 
to that of a non-accreting star of similar spectral type and age. 

The strongest photospheric absorption line that lies in our 
spectral range is Li (6708.1 A). In order to accurately determine 
the depth of the Li line we needed smooth, distinct and clearly visi- 
ble absorption profiles. It is also necessary to have a spectral match 
for each of the accretors in order to estimate the amount of veiling. 
Due to these restrictions, our analysis was limited to four accre- 
tors with distinct Li absorption in their spectra for all 12 epochs, 
IS0126, T45, T30 with spectral types M1.25, M1.25, M2.5 respec- 
tively and T31 with a spectral type of K8. Two non-accretors with 
suitable spectral types and strong Li absorption were used for com- 
parison, T39A (M2) and CHXR20 (K6). We measured the Li EW 
in all epochs for these six objects, and the measurement for three of 
these accretors and their non-accretor spectral match are shown in 
Fig. lAll The non-accretor shows steady deep Li absorption, indica- 
tive of photospheric absorption. All three accretors show shallower 
and more variable Li absorption, which is what we would expect 
when there is variable veiling present in the continuum. (Note: only 
1 1 epochs are considered in the case of ISO 126 due to an increase 
in noise in the second last epoch.) 

In our calculations of the accretion rates we assumed that veil- 
ing was negligible. Since the Li (6708.1 A) absorption line is close 
to the Ha (6562.81 A) emission line, we can consider an estimate of 
the veiling derived at the wavelength of Li to be a good estimate for 
the veiling around the Ha line. We then used this veiling estimate 
to correct our measurements for the Ha EW, as the continuum we 
had measured is actually veiled continuum and not true continuum. 
To find the veiling factor we use the equation 

EWp 

r = - 1 (Al) 

EWv 

where EWp is the photospheric Li EW, and EWv is the veiled 
EW. (For definition of the veiling factor r see Sect. |5). In this case 
the EWp is the EW of the Li absorption line in the reference non- 
accretor, and EWv is the EW of the veiled Li line in the accretor. 



Object 


Log(M) Corrected 


Log(M) Uncorrected 


r 


T45 


-7.42 ± 0.38 


-7.82 ± 0.34 


1.16 


ISO 126 


-8.53 ± 0.36 


-8.65 ± 0.46 


0.32 


T30 


-8.78 ±0.71 


-9.25 ± 0.59 


1.39 


T31 


-6.86 ± 0.36 


-7.06 ± 0.27 


0.47 



We then use the same equation with the calculated veiling factor, to 
derive an unveiled EW for Ha, and used it to recalculate the accre- 
tion rates. This was done for all epoch measurements, and so our 
new accretion rate estimates also take into account the variations in 
veiling. These measurements are presented in Table lAll 

There are a number of sources of uncertainty using this 
method. The veiling factors were derived to give an estimate of how 
much veiling there is in the spectra and hence allow us to evaluate 
the affect of veiling on our results. Ideally to reduce the spread and 
errors in the veiling estimate, a number of lines are used, however 
since we just have a single one available, our calculations can only 
provide a rough estimate of the veiling. The abundance of lithium 
in YSOs decreases with age, so an older star will have shallower 
Li absorption. Previous studies have found stellar clusters, such as 
TW Hydrae and r\ Chamaeleo ntis with ages of ~ 12 Myr, to have a 
range of Li EW of 0.4 - 0.6 A l lMentuch et alj2008h . However at ~ 
2 Myr, Chal is a younger region, where an even smaller range in Li 
EWs can be expected. Similarly, a fast rotator will have a broadened 
Li absorption line, making it appear shallower, and will lead to an 
overestimation of the veiling. In our sample all the accretors have 
narrow Li absorption, so this is unlikely to be a significant factor. 

On average when we take account of veiling within the con- 
tinuum, the accretion rate for all four objects increases by 0.30 dex. 
However more importantly for our results is that the variations in 
the veiling have little effect on the amplitude of variations in the 
accretion rates, there is only a 0.09 dex increase on average in the 
spread. Using this figure we can hypothesize that on average, the 
spread in all the accretors will increase by the same amount. This 
will increase the average spread in accretion rates derived from Ha 
EW to 0.46 dex. The same can most likely be said for the accretion 
rates derived from Call EW, however veiling tends to be less at 
longer wavelengths. Therefore, although the veiling present in the 
continuum may have caused us to underestimate the accretion rates 
and spread slightly, it does not affect our conclusions. 



APPENDIX B: ERROR ESTIMATION 

In order to get an estimate of our measurement errors we varied 
both the sky removal and emission line measurement parameters. 
Specifically we used (i) three different sky continuum estimates, 
the spatial averaged sky continuum and the average sky continium 
± standard deviation (ii) five different window sizes around the cho- 
sen window size given within the text, and (iii) five different win- 
dow positions around the line centre. The ranges in parameters were 
chosen to allow for reasonable measurements of the emission line, 
while still differing from the correct parameters. We varied these 
three sets of parameters were seperately for each measurement of 
the Ha, Ca II and He I emission lines. 

The mean difference in the correct line measurement and the 
perturbed line measurement was found for each of the three param- 
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eter variations. After these mean differences were added in quadra- 
ture, the square root of the sum was taken to get an estimate of the 
total error in each line measurement (see Table [2] and [3} . 

These total error estimates were then carried through to the 
accretion rate estimates, giving the error esimtates in Table[5] 
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